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Abstract. The propagation of 7-rays through metagalactic space is asso- 
ciated with pair creation and subsequent inverse-Compton scattering off 
low-energy background radiation. As a consequence, 7-rays of very high 
energy emitted by remote sources are reprocessed into the window from 
10 MeV to 30 GeV conserving the injected energy. Any cosmologically dis- 
tributed population of 7-ray sources therefore contributes to the diffuse 
7-ray background in this energy band which is well-determined from re- 
cent observations with the Compton Gamma Ray Observatory (CGRO). 
Since the 7-rays trace accelerated particles, the observed flux of diffuse 
7-rays also constrains the global efficiency for particle acceleration. Radio 
galaxies can account for the 7-ray background if their particle acceleration 
efficiency considerably exceeds ~ 18% implying that particle acceleration 
is an essential part of the thermodynamics in these sources. 



1. Introduction 

The microphysics of particle acceleration in the solar system is a fascinat- 
ing realm and provides useful paradigms for particle acceleration in other 
places of the Universe. From the observed flux of cosmic rays it is clear 
that much more powerful particle accelerators than main sequence stars 
and their winds must be at work in our Galaxy and in other galaxies. Su- 
pernova remnants (Ellison et al, 1997) and radio galaxies (Biermann, 1995) 
are among the most interesting putative sources of cosmic rays. The cos- 
mic ray pressure is rather large and compares with that of the interstellar 
magnetic fields implying that cosmic rays are an important ingredient in 
the dynamics of the interstellar medium, and possibly also in the dynamics 
of the intergalactic medium (in superclusters). This imposes strong con- 
straints on the energetics of possible cosmic ray sources from which one 
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can infer the particle acceleration efficiency. For instance, supernova blast 
waves can supply the flux of cosmic rays up to the so-called knee in their 
spectrum at an energy of 10 15 eV only if their particle acceleration efficiency 
is as large as ~ 13% (Drury, 1990) which basically invalidates simple test- 
particle approaches for the description of the acceleration mechanism. This 
has enforced the two-fluid theory for shock acceleration in which the mo- 
mentum flux due to accelerated particles is self-consistently included into 
the dynamics of the shock wave (Achterberg et al, 1984). In this contri- 
bution it is argued that if radio galaxies are the most powerful particle 
accelerators in the Universe responsible for the diffuse isotropic gamma-ray 
background, and possibly also for the ultrahigh-energy cosmic rays up to 
at least 5 x 10 19 eV (the so-called Greisen-Zatsepin-Kuzmin cutoff), they 
require a particle acceleration efficiency even larger than that of supernova 
remnants. 

2. Energy density of the diffuse isotropic 7-ray background 

The energy losses of relativistic particles inevitably lead to 7-radiation by 
which the acceleration sites can be traced. The cumulative flux from all 
unresolved cosmic accelerators appears as an diffuse isotropic flux in an in- 
crtial frame. A diffuse high-latitude 7-ray background has been measured 
with CGRO from MeV to ~ 100 GeV photon energies which very likely 
arises from unresolved extragalactic sources (Sreekumar et al., 1998). The 
resolved sources already contribute ~ 15% of the diffuse isotropic flux and 
show the same average spectral shape as the diffuse background. The back- 
ground spectrum above 30 MeV as determined with EGRET on board 
CGRO is given by 

j AT / F \ - 2 - 1 ±0.03 

=- = (7.32 ± 0.34) x 10- 9 cm' 2 s" 1 MeV" 1 sr" 1 . 

dE K ' V 451 MeV/ 

(1) 

The spectrum continues into the MeV range. Most of the background flux 
below 10 MeV can be explained by the contribution of Supernovae la from 
the era of galaxy formation at Zf = 2—5 (The et al., 1993). The uncertainties 
in this thermal component of the diffuse 7-ray background are not very 
large, the known element abundances tightly constrain the total number of 
SNIa, and hence the total energy flux in 7-rays, in the Universe leaving little 
room for other sources contributing in the MeV range. The uncertainty in 
the cosmic star formation history due to possible dust enshrouding of the 
earliest galaxies should amount to less than a factor of three. Integrating 
the above non-thermal background spectrum between 10 MeV and 30 GeV 
(the flux at higher energies must be considered somewhat unreliable due to 
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the strongly decreasing aperture of EGRET) one obtains the energy density 
u 7 = (5.03 ± 0.64) x 10~ 6 eVcm- 3 (2) 

or 

ft 7 /i 2 = (2.09 ± 0.27) x 10~ 10 (3) 
where h = H o /W0 km s^ 1 Mpc -1 parametrizes the Hubble constant. 

3. Metagalactic gamma-ray absorption and production 

Calorimetry of the particles accelerated throughout the Universe requires 
to know whether all 7-rays emitted by a source also reach the present-day 
observer. However, 7-rays of energy E can interact with low-energy photons 
of energy e from the diffuse isotropic background over cosmological distance 
scales I producing electron-positron pairs 7 + 7 — > e + + e~ , if their energy 
exceeds the threshold energy 

2(m c c 2 ) 2 fl + z\~ 2 { E \- 1 , s /N 

where denotes the cosine of the scattering angle (Gould & Shreder, 1966; 
Coppi & Aharonian, 1997). The 7-ray attenuation e~ T due to pair produc- 
tion becomes important if the mean free path A becomes smaller than I, 
i.e. if the optical depth obeys r = l/X > 1. For the computation of r one 
first needs to know the pair production cross section 



2/3(/3 2 -2) + (3-/3 4 )ln(i±|) 



(5) 



where j3 = y/1 — I/7 2 with 7 2 = e/e t h, and where <tt denotes the Thomson 
cross section (Jauch & Rohrlich, 1976). Then one needs the geodesic radial 
displacement function dl/dz = 77- [(1 + z)E(z)]~ 1 to compute the line in- 
tegral from z = to some z = z . For a cosmological model with = 1 
and A = the function E(z) simplifies to (1 + z) 3 / 2 . Hence one obtains the 
optical depth 

rz dl r+1 1 — fl f°° 
Tjj(E, Zo) = / dz— / d[i — - — / den h (e)(l+z) cr 77 (^,e, fi,z) 

JO dZ J-l 2 J 6th 

c f z ° r 2 x f°° 

= — / dz(l + z) 1/2 dx- den h (e)a^(E,e,x - l,z) (6) 
Ho JO JO 2 J e . 



adopting a non-evolving present-day background density nt>, i.e. n' b (z, e')de' = 
(1 + z) 3 n\ ) (e)de where the dash indicates comoving-frame quantities. The 
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Figure 1. The diffuse isotropic microwave-to-ultraviolet background: the solid curve 
shows a 10th order polynomial interpolation of observational data ((Fixsen el al, 1998; 
Mannheim, 1998; Franceschini el al., 1998), and references in (Madau & Pinney, 1996)). 

simplifying assumption that the photon density transforms geometrically 
corresponds to a situation in which an initial short burst of star formation 
at Zf > z produced most of the diffuse infrared-to-ultraviolet background 
radiation. Fig.l shows the spectrum of the low-energy diffuse background 
used to solve Eq.(6) numerically. 

Figure 2 shows the resulting t(E, z) = 1 (omitting the subscript here- 
after) curve for the microwave-to-ultraviolet diffuse background spectrum 
shown in Fig.l. It is obvious that 7-rays above ~ 10 — 50 GeV can not 
reach us from beyond redshifts of z = Zf = 2 — 4. Higher energy 7-rays can 
reach us only from sources at lower redshifts (e.g. 7-rays with energies up 
to 10 TeV have been observed from Mrk 501 at z = 0.033 in accord with 
Fig.2 (Mannheim, 1998)). 

Corollary I: If the DIGB originates from unresolved sources distributed 
in redshift similar to galaxies, its spectrum must steepen above ~ 30 GeV 
due to j-ray pair attenuation. 

Here is has been tacitly assumed that the 7-rays which have turned 
into electron-positron pairs do not show up again. This is, in fact, not 
quite true, since the pairs are subject to inverse-Compton scattering off the 
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Figure 2. The 7-ray horizon t(E,z) = 1 for the low-energy background spectrum 
shown in Fig.l. Cosmological parameters are h = 0.6, Q = 1, and SIa = 0. For a general 
discussion of pair attenuation, see reference (Biller, 1996). 

microwave background thereby replenishing 7-rays. The 2.7 K background 
is more important as a target than the shorter wavelength background, 
since there is no threshold condition for Thomson scattering contrary to 
pair production and since 2.7 K photons greatly outnumber the latter. The 
inverse-Compton scattered microwave photons turn into 7-rays of energy 



conserving the energy of the absorbed 7-ray which corresponds to a constant 
E 2 dN/dE, i.e. the expected slope of the differential spectrum is about -2 
(-2.1 observed). A negligible amount of energy is lost to lower frequency 
synchrotron emission, if magnetic fields are present in the interagalactic 
medium. 

Corollary II: Energy conservation in the reprocessing of "y-rays from 
higher to lower energies by pair production and subsequent inverse-Compton 
scattering produces an approximate dN/dE oc E~ 2 power law DIGB be- 
tween ~ 10 MeV and ~ 30 GeV. 
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Figure 3. Extragalactic background radiation. The solid curve interpolates observational 
data. Whereas AGN dominate in the X-ray regime, their longer wavelength emission 
(dotted line from (Sanders et al, 1989)) is surpassed by the stellar and dust emission 
from early galaxies. 



4. Particle acceleration efficiency 

The energy density of the accelerated particles in the sources (u acc ) must 
obey u acc > u 7 for consistency. Defining the radiative efficiency 



£rad = — = Mm 



it 



(8) 



where t a and t c denote the adiabatic and energy loss time scales, respec- 
tively, the energy requirement for the sources is minimized for £ ra d = 1. 
In extragalactic radio sources, £ ra d ~ 1 corresponds to very large electron 
(7e ~ 10 5 ) or proton Lorentz factors (7 P ~ 10 10 ) consistent with the non- 
thermal 7-ray spectra extending above 10 GeV. 

The resolved extragalactic 7-ray sources belong to the so-called blazar 
subclass of radio-loud active galactic nuclei (AGN) in which one observes 
jets emerging from accreting supermassive black holes preferentially at 
small angles to their relativistic velocity vector. If the entire class of radio- 
loud AGN is responsible for the DIGB, the particle acceleration efficiency 
is constrained by 

Ucc = ^£ = i (9 ) 
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where Uj denotes the total (kinetic + magnetic + randomized relativistic 
particle) energy density in extragalactic jets. Since u 7 is known (Eq.(l)), 
a determination of Uj is necessary to further constrain This energy 
density can be inferred from the relative strengths of the various bumps in 
the overall diffuse background spectrum (Fig.3) in the following way. 



5. Origins of the diffuse isotropic background radiation 

The microwave bump is the well-known signature of the big bang at the 
time of decoupling with its energy density given by the Stefan-Boltzmann 
law = ctT 4 . The bump in the far- infrared is due to star formation in 
early galaxies (at z = Z{), since part of the stellar light, which is visible 
as the bump at visible wavelengths, is reprocessed by dust obscuring the 
star-forming regions. The energy density of the two bumps can be inferred 
from the present-day heavy element abundances. Heavy elements have a 
mass fraction Z = 0.03 and were produced in early bursts of star formation 
by nucleosynthesis with radiative efficiency e = 0.007 yielding 

UM ~ TT _~6xlO 3|__j^__j eV cm3 (10) 

About half of the energy is contained in either bump. It has been shown 
recently by a number of groups that probably all galaxies (except dwarfs) 
contain supermassive black holes in their centers which are actively accret- 
ing over a fraction of i a gnA* ~ 10~ 2 of their lifetime implying that the 
electromagnetic radiation released by the accreting black holes amounts to 

e accr-^bh ^agn -. . 1n -4 T —3 

«accr T^M 1 ~ Cm ' ' 

adopting the accretion efficiency e accr = 0.1 and the black hole mass frac- 
tion M bh /M* = 0.005 (Rees & Silk, 1998). Most of the accretion power 
emerges in the ultraviolet where the diffuse background is unobservable ow- 
ing to photoelectric absorption by the neutral component of the interstellar 
medium. However, a fraction of u x /ubh ~ 20% (from the average quasar 
spectral energy distribution (Sanders et al, 1989)) shows up in hard X- 
rays due to coronal emission from the accretion disk to produce the diffuse 
isotropic X-ray background bump with u x - 2.8 x 10~ 5 eV cm (Gruber, 
1992). 

Jets with non-thermal (7-ray) emission show up only in the radio-loud 
fraction £ r i ~ 20% of all AGN and their kinetic power roughly equals the 
accretion power (Rawlings & Saunders, 1991). Hence one obtains for the 
energy density in extragalactic jets 

£rl \ ( £rl 



0.2 



(Jj^J 2.8 x 10~ 5 eV cm" 3 (12) 
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Substituting this value into Eq.(9) one obtains a limit for the acceleration 
efficiency 

^°- i8 (tjHi§r <i3) 

which has to be compared with the 13% efficiency required for supernova 
remnants. Since a radiative efficiency much lower than unity must be con- 
sidered realistically (e.g., to account for the pdV work of the jets against a 
surrounding intracluster medium), the acceleration efficiency may have to 
be as large as 30% or more. Possible escapes to this conclusion are (i) the 
radio-loud fraction was much larger at high redshifts than at low redshifts, 
(ii) there is another contributor to the gamma-ray background, such as the 
decay of topological defects (Sigl, 1996), or (hi) the flux of the DIGB is 
actually lower due to the foreground from a galactic halo (Dixson et al., 
1998). 

6. A remark on cosmic rays from radio galaxies 

An extragalactic cosmic ray component seems to be evident from the ob- 
served change in slope of the local spectrum above 3 x 10 18 eV. The slope of 
this extragalactic component is much steeper than E~ 2 consistent with the 
steepening due to energy losses for a cosmologically distributed and evolv- 
ing source population, such as extragalactic radio sources (Rachen et al., 
1993). The total energy density in the loss-steepened extragalactic cosmic 
ray component and the DIGB are of the same order of magnitude, as re- 
quired if the 7-rays trace neutral pions and the cosmic rays trace neutrons 
produced in photo-meson production events. 

7. Summary 

The particle acceleration efficiency of the jets in radio galaxies must be 
larger than ~ 18% if the 7-rays from these jets are to explain the DIGB. 
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